Matter Coupling to Strong Electromagnetic Fields in Two-Level Quantum Systems 

with Broken Inversion Symmetry 



ON 

o 
o 

(N 

c 

a 

(N 



O.V. Kibis 1 ^ G.Ya. Slepyan 2 , S.A. Maksimenko 2 , and A. Hoffmann 3 

1 Department of Applied and Theoretical Physics, Novosibirsk State Technical University, 
Karl Marx Avenue 20, 630092 Novosibirsk, Russia 
^Institute for Nuclear Problems, Belarus State University, Bobruyskaya St. 11, 220050 Minsk, Belarus 
3 Institut fur Festkorperphysik, Technische Universitdt Berlin, Hardenbergstrajie 36, D- 10623 Berlin, Germany 

We demonstrate theoretically the parametric oscillator behavior of a two-level quantum system 
with broken inversion symmetry exposed to a strong electromagnetic field. A multitude of reso- 
nance frequencies and additional harmonics in the scattered light spectrum as well as altered Rabi 
frequency are predicted to be inherent to such systems. In particular, dipole radiation at the Rabi 
frequency appears to be possible. Since the Rabi frequency is controlled by the strength of coupling 
electromagnetic field, the effect can serve for the frequency-tuned parametric amplification and gen- 
eration of electromagnetic waves. Manifestation of the effect is discussed for Ill-nitride quantum 
dots with strong build-in electric field breaking the inversion symmetry. Terahertz emission from 
arrays of such quantum dots is shown to be experimentally observable. 
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The resonant interaction of quantum systems with 
strong electromagnetic field [J [J is permanently in the 
focus of interest, both due to high methodological value 
of arising problems and their direct relation to current 
appliedprojects, such as new generation of high-efficient 
lasers Q, laser cooling of atoms Q, development of ba- 
sis for quantum information processing 0], etc. One of 
bright manifestations of the strong field-matter coupling 
is the Rabi effect Q: oscillations of the level population 
in a quantum system exposed to a monochromatic elec- 
tromagnetic wave. The simplest physical model leading 
to harmonic Rabi oscillations is a two-level symmetrical 
quantum system placed in a given classical single-mode 
electromagnetic field Incorporation into this sim- 
plest model additional physical factors results in many 
nontrivial effects. For example, accounting for the quan- 
tum nature of light leads to the concept of radiation- 
dressed atoms [l[ and the 'collapse-revivals' phenomenon 
in the population dynamics of a system exposed to coher- 
ent light. Time-domain modulation of the field-matter 
coupling constant [f| , local- field effects in nanostructures 
0, 0] and phonon- induced dephasing [§] provide new pos- 
sibilities for the control of the Rabi oscillations. Many 
interesting effects manifest themselves in more complex 
systems, such as two coupled Rabi oscillators [§], Rabi os- 
cillators based on superconducting electrical circuits 
and systems where Rabi oscillations are strongly influ- 
enced by intraband motion of quasi-particles [lfj . In the 
given Letter we investigate the role in the Rabi effect of 
a new physical factor — violation of the inversion sym- 
metry. Inherent in many quantum systems, the factor is 
ignored by the conventional physical model 0] of Rabi 
oscillations. Meanwhile, our theoretical analysis predicts 
pronounced manifestation of the violation in a set of ob- 
servable effects in different physical systems. 

Let us consider a two-level quantum system with \a) 



and 1 6) as excited and ground states, respectively. Let 
the system interacts with a classical linearly polarized 
monochromatic electromagnetic field. The system is de- 
scribed by the Hamiltonian H = Hq + Hint- The free- 
particle Hamiltonian, written in the basis of these two 
states, is TCq — Hu>oa z /2, where a z is the Pauli matrix 
and wo is the resonant frequency of the two-level sys- 
tem. The interaction Hamiltonian is expressed in terms 
of the amplitude E and frequency lo of the driving field 
by Hint = — Edcos(ojt), where d is the electric dipole mo- 
ment operator and d^ = (£|d|j) are its matrix elements. 
The critical assumption, which distinguishes the systems 
being considered from standard ones, is the violation of 
the inversion symmetry. Since in that case the states \a) 
and 1 6) do not possess a certain spatial parity, the diago- 
nal matrix elements of the dipole moment operator prove 
to be nonequivalent, d QQ ^ dbt, dictating thus physical 
effects described hereafter. 

We shall seek the solution of the Schrodinger equa- 
tion with the Hamiltonian described above in the form 
of \vb) = C a (t)\a) + Cb(t)\b). Substituting this expression 
into the Schrodinger equation, we arrive at the equations 



ihC a = 
iHCb = 



HuJn 



flWn 



- Ed aa cos(arf) C a - Ed ah cos(wt)C 6 , (1) 
- Edf, b cos(uit) Cb - Ed ba cos(wt)C a .(2) 



In order to solve these equations with respect to C a 
and Cb, we first rewrite them for modified amplitudes 
c a ,b = C ai bexp[±iL0 t/2 - i<j> atb (tjj,t)], where (pj(ui,t) = 
Edjj sin(ujt) / hco , and signs ± correspond to indexes a 
and b, respectively. The system of equations (H])-© is 
then reduced to the form as follows 



hc a = iE* ff d ab cos(^)c fc e Mot , 
hc b = iE c ffd fea cos(^)c a e~ lli ' ot 



(3) 
(4) 
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where the effective electric field strength is 

oo 

E cff (i) = Ee - lKsin ^ = E ]T J »(«) e 



} inujt 



(5) 



J n (/c) is the Bessel function of the first kind, and 
K = E(df,f, — d aa )/huj is the symmetry violation pa- 
rameter. Let us stress that Eqs. I©-® are analogous 
to standard equations of two-level system Q with the 
only difference that the driving field amplitude E is re- 
placed by the effective amplitude (JSJ). That the effec- 
tive amplitude E e g(t) is time-dependent causes the sys- 
tem to become a parametric oscillator with full set of 
intrinsic properties, unusual for standard Rabi oscilla- 
tors. In particular, a multitude of resonant frequencies 
lu = luq/h, n = 1,2,3,... appears in the system. Fur- 
ther we shall assume that the driving field frequency u> is 
in the vicinity of the frequency loq /m (m-th resonance) 
and the conditions \^id ao mJ m (K) / kH(ujq — mu>)\ > 1 and 
|Ed a |,n J n (n) / K,h(ujQ — nuj)\ < 1 (n ^ m) are fulfilled, 
allowing us to neglect interaction of the two-level system 
with all harmonics ([5]) other than harmonics mil. Then 
Eqs. (J3j) — (J4j) resume the form analogous to equations for 
symmetric two-level systems and can easily be solved by 
invoking the rotating- wave approximation Q . Assuming 
the nondiagonal dipole matrix elements to be real- valued, 
we obtain 



C a {t) = |C Q (0) 

.Sl R 



cos 



Sit 

T 

(Sit 



iA /Sit 
— sin — 
SI V 2 



+t-^-C b (0) sin(y) I e -w/2 e i* a (*,t) ( (6) 



C b (t) = C 6 (0) 



+z^C a (0) sin(^) \ e ^t/2^,t) ^ (?) 



(Slt\ iA (Sit 

C0S ItJ + iiHy 



(Sit 



where the parameters A 
and the Rabi frequency 



V 2 



muj, SI 



Sl R = 2Ed ab mJ m (n) / nh 



(8) 



are written in the vicinity of m-th resonance. Since we 
consider the strong coupling regime, the driving field E 
is assumed to be sufficiently strong, SIrt 3> 1, to ne- 
glect impact of the linewidth H/t. Note that the Rabi 
frequency in systems with broken inversion symmetry 
(jHJ) shows nonmonotonic dependence on the driving field 
strength and even can turn zero at those values of E 
which correspond to roots of the Bessel function J m (re). 

It is well known that in symmetrical two-level systems 
the Rabi effect manifests itself in the power spectrum 
of the scattered light by peaks centered at the incident 
light frequency u and at the displaced frequencies u ± SI 
(Mollow triplet, [HI]). I n order to reveal peculiarities of 
the scattered light induced by the symmetry violation, 



it is suffice to consider the process in the framework of 
classical electrodynamics. It allows analyzing electronic 
subsystem, in respect to irradiation of electromagnetic 
waves, as a classical dipole with the oscillating dipole 
moment d(i) = (ip\d\ip) . For definiteness, we identify 
initial time as that corresponding to system being in the 
excited state, i.e., C o (0) = 1 and C&(0) = 0. Then, 
substituting wave function with coefficients ©-(O 
into expression for the dipole moment and omitting the 
time-independent terms, we arrive at the expression 



d(t) = (d aa - d 66 ) 



SI 



R ASH 



x S J m _ n (ft)e 
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(9) 



for the dipole moment in the vicinity of m-th resonance. 
As follows from Eq. ©, the radiation spectrum consists 
of a singlet at the frequency SI, and an infinite sequence 
of triplets with the frequencies no;, nuj±Sl (see Fig. 1). It 
should be emphasized that the frequency multiplication 
is the general property of systems with broken inversion 
symmetry and can be used for the high harmonic genera- 
tion [l2| . Amplitudes of harmonics of the triplets rapidly 
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FIG. 1: Schematics of the light scattering in a two-level sys- 
tem with broken inversion symmetry in the vicinity of one of 
the possible resonances u = Uo/m. 

decrease with increasing norm, while the singlet SI am- 
plitude depends on neither n nor m. Radiation of the 
dipole © in the vicinity of each resonance is character- 
ized by its own Rabi frequency ([5]) decreasing with the 
resonance number m increase. In systems with inversion 
symmetry diagonal elements of the dipole moment are 
identical. As a result, higher resonances (m > 1), higher 
triplets (n > 1), and the singlet vanish. In that case, 
coefficients ©-([7]) and expression for the dipole moment 
([9]) coincide with the solution of the problem presented 
in Ref. @. 

Further we confine consideration to systems with a 
weak violation of the inversion symmetry, when the con- 
dition |k| -c 1 holds true. In such systems the asymmetry 
effect is expected to be most pronounced in the vicinity 



3 



of the first resonance m — 1 (ui ~ ljq). If then we restrict 
ourselves to the most interesting case of resonant driving 
field, when A = and ft = Q R , we reduce © to the 
expression 



d(t) 



cos(Q R t)-d ab sm(uj t)sin(tt R t) , (10) 



where Rabi frequency (JSj> takes its conventional form 
fin = Edab/h. In typical quantum systems the Rabi 
frequency is much less then the driving field frequency 
lu = luq. Therefore, in addition to the high-frequency 
harmonics loq ± fl R (second term in the right-hand part 
of Eq. (HHJ)), violation of the inversion symmetry leads to 
the irradiation of low- frequency electromagnetic waves at 
Rabi frequency (first term in the right-hand part of Eq. 
(|10p ) with the time- averaged radiation intensity 
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The quantity |d aa — d b| is further referred to as effective 
dipole moment. Note that the frequency of the scattered 
radiation £l R depends only on the driving field E and 
does not depend on both the frequency of the scattering 
system loq and the frequency of incident light lu. 

The simplest quantum system devoid of the inversion 
center is a Rydberg hydrogen atom imposed to a homo- 
geneous static electric field £ . Assuming the atom to be 
in ground state, we find that the Rabi frequency is ex- 
pressed in terms of the Bohr radius as = h 2 /m e e 2 by 
Q R = 4(2/3) easE/h. Correspondingly, for the elec- 
tric field £ much less than the intra-atomic electric field 
e/a%, the intensity of the radiation at the Rabi frequency 
is determined by Eq. (fTTj) with the effective dipole mo- 
ment given by |d Qa - d b6 | = (l/8)(4/3) n a%£. This ex- 
pression is applicable not only to hydrogen atom but also 
can be used for estimating parameters of the emission at 
the Rabi frequency in arbitrary quantum systems with 
broken inversion symmetry. For that aim, the Bohr ra- 
dius as should be replaced by the characteristic linear 
extension of the system. Therefore, the radiation inten- 
sity (JTTJ) rises with the increase of the system size. In that 
connection, confined semiconductor nanostructures with 
discrete energy spectrum and linear extension multiply 
exceeding the atomic size, quantum dots (QDs) 0, can 
serve as prospective systems for the effect observation. 

In quantum dots, driving electromagnetic field trans- 
fers electrons from valence band into conduction band; 
what is why in QDs ground state \b) corresponds to 
the absence of free carriers while first excited state \a) 
is the state with electron in the conduction band and 
hole in the valence band. In that case, Hluo is approx- 
imately equal to the semiconductor bandgap, and the 
Rabi frequency is determined by the standard expression 
Q R = Ed c „/7i with the dipole matrix element d cv cor- 
responding to interband transitions. Among a variety 



of different types of QDs, nitride-based confined struc- 
tures seem to be most promising for the effect observa- 
tion (for parameters of Ill-nitrides A1N, GaN, and InN 
see Ref. [13J). Indeed, GaN and similar Ill-group nitride 
alloys have a hexagonal (wurtzite) structure. As a con- 
sequence of giant piezoelectric effect inherent in hexag- 
onal crystals, QDs based on structures AIN/GaN and 
InN/GaN have a strong built-in strain- induced electric 
field with strength £ of several MV/cm [13, [Hj]. Due to 
the strong electric field, conduction-band electrons and 
valence-band holes in the QDs get spatially separated 
[HI, 12 1 1 an d the effective dipole moment of a Ill-nitride 
QD is given by simple relation |d aQ — d^ ~ el, where I is 
the QD height. Then, for typical I of several nanometers, 
the effective dipole moment of nitride QDs is estimated 
as |d aa — dbb\ ~ lODebye, what is tens thousands times 
as large as the effective dipole moment of hydrogen atom 
in the same electric field. 

Since the Rabi frequency depends on the driving field 
strength E, the broken inversion symmetry-induced low- 
frequency singlet fiij in the dipole emission spectrum 
can be used for the tunable generation of electromag- 
netic waves at Rabi frequency with intensity (fTTj) . This 
problem is especially challenging for frequency ranges 
where traditional methods either fail or inefficient, such 
as terahertz domain. Since this domain lies between ra- 
dio and optical frequency ranges, neither optical nor mi- 
crowave techniques are directly applicable for generating 
THz waves. Therefore, a search for effective THz ra- 
diation sources is one of most excited fields of modern 
applied physics 1|| 1^, 2(J. One of the latest trends to 
fill the THz gap is using nanostructures as THz emitters 
and detectors. The quantum cascade THz transitions in 
QD systems 2l|, |22j and different electron mechanisms 
of THz emission from carbon nanotubes 23j, |2J, |25[ have 
been proposed and are being actively studied. Thus, the 
proposed mechanism of THz emission fits well the current 
tendencies in nanophotonics. 

Parameters of the THz emission from Ill-nitride QDs 
can be evaluated in the following way. Using the esti- 
mate dgj, ~ lODebye for the QD interband dipole mo- 
ment [13j . we obtain that the Rabi frequency turns out 
to be lying in the THz range at relatively weak strength 
of the driving field, E ~ 10 5 V/cm. Then, substitut- 
ing the previously obtained value of the effective dipole 
moment d aa — d^l ~ lODebye into Eq. (JTTJ) , for the 
typical QD area ~ 10~ 12 cm 2 we arrive at the radia- 
tive intensity per single QD ~ 10 -11 W/cm 2 in the THz 
range. This estimate multiples as N 2 for array of N 
identical QDs with the lateral extension less than the 
Rabi wavelength = 2nc/Q R . In such an array, all 
QDs emit waves in phase. The THz-range wavelength 
A_r ~ 10 -2 cm restricts the array lateral size. Taking 
into account that the typical density of nitride QDs is 
~ 10 11 cm -2 , we can estimate the number of QDs in such 
an array by N ~ 10 7 — 10 s . Therefore, the THz emission 
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power from the submillimeter-sized QD array may ap- 
proach the micro- Watt level, which is characteristic for 
the resonant-tunneling diodes based on carbon nanotubes 
(2^ . Certainly, state-of-the-art THz quantum cascade 
lasers demonstrate essentially larger output power [27 1. 



However, array of QDs with broken inversion symmetry 
provides smooth frequency tuning by the variation of the 
driving field intensity. The reducing impact of inhomo- 
geneous broadening is beyond of the estimate and will 
be considered elsewhere. As to the homogeneous mecha- 
nisms of the line broadening, like phonon scattering, they 
are accounted for phenomenologically in the linewidth 
H/t and do not impose any specific restrictions on the 
observability of the effect, but just influence the criterion 
of the strong light-matter coupling regime. Thus, the 
upper estimate presented of the radiation output allows 
proposing the broken inversion symmetry-induced mech- 
anism of the radiation for the development of novel-type 
THz emitters based on QDs avoid inversion symmetry. 
Obviously, parametric amplification of THz radiation in 
QD arrays is also possible and can be applied for THz 
detecting. 

Once more type of asymmetrical two-level quantum 
systems where the discussed effect can be observed is su- 
perconducting quantum qubits @, HI]. Away from the 
charge degeneracy point of the superconducting quan- 
tum circuits [1| formed by Josephson qubits in microstrip 
resonators, the qubit symmetry is broken and the Hamil- 
tonian of the system (see Eq. (16) in Ref. 0]) takes the 
form analogous to that used in our analysis. At that, the 
resonant transition frequency amounts to ~ 10 GHz while 
the Rabi frequency lies in the range ~ 100 MHz. Inten- 
sity of low-frequency line in the spectrum of Rabi oscil- 
lations can be controlled by the changing of dc gate volt- 
age. Chiral nanostructures, including chiral carbon nan- 
otubes [1^], should also be noted as prospective systems 
for observation of the effect, since the chirality breaks 
the inversion symmetry. New interesting effects are ex- 
pected for asymmetric two-level oscillators placed inside 
band-gap structures like photon crystal, microcavity or 
nanoantenna. For example, by corresponding manipula- 
tion of the photonic crystal parameters, the frequencies 
lj and loq can be chosen laying in the band gap while the 
frequency fl lies in the transparency band. In that case, 
the strong coupling regime is realized for the pump field, 
while at the frequency £1 the photonic crystal serves as 
an antenna transforming near field into far field. Simi- 
lar effect appears in microcavity with resonant frequency 
close to u> and leaking modes at the frequency Q. 

In the paper we have assumed that the incident light 
pulse is much longer than periods of both incident and 
scattered electromagnetic field. When a collection of os- 
cillators with broken inversion symmetry is illuminated 
by an extremely short pulse, the difference between am- 
plitudes E and E e $ in Maxwell-Bloch equations may re- 
sult in the failure of the area theorem and in the related 



effect of the carrier- wave Rabi flopping. Such an effect 
is observed in semiconductors when the Rabi frequency 
becomes comparable to the band gap frequency [301 ] . 

To conclude, let us stress that violation of the inversion 
symmetry is common property of quantum oscillators of 
different physical origination and, consequently, the ef- 
fect predicted is expected to manifest itself in spectral 
characteristics of different optical processes, e.g. reso- 
nant fluorescence of molecules. 
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